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ABSTRACT 


An algorithm has been developed to predict the 
thermal behaviour of a salt gradient solar pond . The 
developed model incorporates a detailed representation 
of the losses from the surface and the bottom of the pori 
and uses hourly metrological data for Delhi to predict the 
pond performance. The effects of x-zater clarity, consecutive 
days of cloudiness, evaporation loss and different rates of 
heat removal on a solar pond's performance have also been 
studied. Validity of assuming a constant transmission 
coefficient for a solar pond throughout the year has also 
been verified . It is concluded that though water clarity 
greatly affects the overall performance of a solar pond, 
variable transmission coefficients have very’’ little influence 
on it. The evaporation loss from the surface of a solar pond 
which is quite significant during the months of May’-, June £: 
July every year, does not allov.' high LCZ temiperat-ures to be 
reached. Stud-y of several consecutive cloudy day’-s sho'ws that 
solar pond still offers a relatively reliable solar collector 


cum-storage system 
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CHAPTER 1 


INTRODUCTION 


1 .1 ENERGY; CRISIS AND RENEWABLE SOURCES 

It is well knov;n that energy represents the source 
of activities of himan beings and nature. It was the discover 
of new sources of energy’ like coal, petrole'um products and 
hydroelectricity which provided a tremendous impetus to the 
industrialisation and it was the quantitative availability of 
these energq^ resources v/hich mainly contributed to the pace 
of development of the, nov.?, industrialised nations. But, 
unfortunately, the pattern of consumption in the first half 
of this century has been, rather, wasteful as no attempt seems 
to have been made to conserve energy with a foresight into the 
future , 

The oil embargo of the early seventies shifted the 
focus to the renewable sources of energy,’-. There is hardly any 
disagreement on the importance of harnessing renewable sources 
of energp', especially solar energy in the immediate future. 

At a time when world is becoming more and more, aware of the 
problems of increased pollution and mounting cost, this clean 
source of energy deserves a serious attention in the world 
energy’- plans . 
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1.2 SOLAR ENERGY 

Sun [l/2] the source of solar energy, is a sphere of 

intensly hot gaseous matter. The rate of energy generation, 

12 

due to fusion reactions, on the sun is around 12x10 Q 

(10 = lO'^^ Joules) . Out of this only a very small portion, 

around 5300 Q enters fne earth's atmosphere. Even this small 

fraction is very huge compared to the total v/orld energy 

consumption of 0.3 Q per annum. However all this energy- that 

enters the atmosphere, does not reach the surface of the earth. 

VJhile at a point just outside the atmosphere, the energy' 

2 

intensity is 1353 w/m , the average amount received on the 

2 

earth is just about 690 Iv/m . This value may vary ivu.th locations, 
seasons i.e., the day of the year, end the time of the day. 

It is observed that the sites lying between laoitucle 
of 30 ®N and 30 ®S have sufficiently high solar intensities even 
during winter. India, with its extreme latitude of S°N and 
32 ®N is fortunately very^ well placed for solar energy 
utilization. Actual solar radiation measurement in various 
ESCAF (Economic and Social Commission for Asia and uhe Pacific) 
countries have been used to produce charts of total solar 
radiation per unit area on unit reception basis. The attached 
chart of Fig. (1.1) shov/s that India receives the ma>umuiTi solar 
radiation. Average reception for India works out to be 2100 
kv'Jh per year per square meter area. 
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Main problems being faced by the solar energy 
utilization programme are its low intensities t'nerby requiring 
large collector areas, lovj overall system efficiencies, and 
the need of auxiliary equipments due to its intermittent 
availability . 

1.3 CONCEPT BEHIND A SOLAR POND 

The separation of the collecting and storage units 
in the conventional solar system permits the efficient design 
for each unit. Although the efficiency of such systems is 
quite high so is the unit cost. These systems can not be 
scaled up for long term, large energ*f storage. The large 
collecting areas ana large storage units will be prohibitively 
expensive. The long term energ^^ storage calls for a different: 
approach. For solar energg^ to be competitive, the system must 
be designed where large collecting areas or large storage 
capacity or both should be cheaply available . Solar pond is 
one of the concepts where both large collecting areas and 
storage capacity are in expensive. 

1 .4 WHAT IS A SOLAR POND ? 

Solar ponds [3,4J are proven, viable, large area 
collectors of solar energg^. These innovative systems collect 
and store solar energy simultaneously. In a water tank with a 
black bottom., solar insolations are absorbed at the bottomi and 
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the water in contact is heated up. Heat is transferred 
from lower layers to upper layers by convection of hot v/ater 
which being lighter than the cold v/ater, shifts to the top. 
Energy is lost to the ambient by the hot upper layer, unless 
means are devised to stop this, most of the energy will be 
■wasted. Thus to store energy^ either glazing should be used 
or convection must be stoped, i.e., "ways must be devised to 
stop this upv;ard movement of "water. 

The unique characteristic of salinit-y gradient pond 
is that, for a large part of their depth, they are non- 
convective . It is just this stability against transfer of 
heat by thermal convection that make possible effective 
collection and storage of solar energy?. Stability against 
convection in saliriity gradient ponds is the result" of a down 
■ward increase of salt concentration, from close to zero at -the 
surface to typically 20-25% near the bottom. The variation in 
salinity makes the "water denser at the bottom even when it is 
heated to a relatively high temperature by the transmitted 
solar radiation. Being denser, it remains at the bottom and 
so its heat is not lost by convection into the surface layer 
and utimately to the atmosphere. Salinity variation inside 
the pond has little effect upon either the solar energy? 
collection or the storage but it stabilizes the pond against 
convection and thus against heat loss produced by convection 
in ordinary ponds and lakes . 
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Given this stability, the only natural process by 
which heat can be transferred upward from the bottom of a 
salinity gradient pond is heat conduction, water being opaque 
to low temperature radiations. Kon-convective water thus 
serves as a partially’- transparent thermal insulator. One 
meter of non-convective ’.va-cer has about the same insulating 
value as six cms of foam, or about five times the insulating 
value of double glass as commonly’- used for V7inacw and solar 
energy/ collectors. 

The fact that water is not entirely transparent to 
solar radiation has important consequences for the pond 
performance . The infrared and red components of the 
spectrum are absorbed near the surface and a little over one 
third of the incident solar energy’ penetrates to a depth of 
one meter in clean water (fresh or saline) . This fact limits 
the collection efficiency’ of solar ponds. However, as 
compared with glass insulated solar collectors, they have the 
com.pensating advantage that there is no minimum threshold 
radiation level for collection end they collect the same fraction 
of the incident energy’ ho'wever low the solar input rna-;/ be . 

1.5 CLASSIFICATION OF SOLAR PONDS 

The word solar pond in its generalised sense is used 
to mean a sy’stem where v/ater is contained in a man m.ade pond 
or natural cavity inside earth. The water simultaneously’ 
collects and stores energy. The water does not flow past an 
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absorbing surface as is the case of conventional water 
heaters. Solar ponds are classified as follows [S]. 

(1) SHALLOW SOLAR POND (SSP) 

It is a large water tank containing water about 
5-10 cm in depth. It has a black bottom and a glasin-g 
at the top to reduce convective losses . The glazing 
is in touch v/ith the v.?ater . The sides of the tank 
are v^eli insulated. Most of the solar insulation is 
absorbed by the blackened bottom and only very little 
by the successive layers of water. The water is 
heated by convection. The pond is filled in the 
morning by cold water and the hot water is drained out 
in the evening. The hot v:ater is either used or stored 
in a well insulated tank and is returned to the pond 
next morning to collect more energy. SSP provides hot 
w’ater in the range of 40 to 60 ®C, There is no long 
term energy storage, it can at best provide short term 
storage . 

In some of the designs, the v;ater is filled into 
a plastic bag which rests on a well insulated black 
surface . 

(2) SALT GRADIENT SOLAR PONDS (SGSP) OR NON -CONVECTIVE 
SOLAR PONDS OR SOLAR PONDS (SP ) 

This is the concept usually known as the solar ponds . 
Here water is contained in a big open tank, artificial 
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or natural, 2 tc 3 m deep. The vjater contains salt, 
mostly NaCl or Mg CI 2 such that the concentration of 
the dissolved salt increases with depth . The solar 
insolation is absorbed by v;ater at different layers 
but mostly by the dark coloured bottom of the ponds . 

A critical value of the salt density gradient is 
needed tc off-set the thermal density gradient set 
up by the hot v.’ater in the bottom layers. Once the 
concentration of dissolved salt is greater than the 
critical value, convection will stop and the hot 
water stays at the bottom. The only mode of 
transmission of neat tc upper layers is conduction, 
which is a much slower process. The water thus acts 
as insulation and heat remains trapped in the bottom 
layers. Under good operating conditions, the storage 
layers can reach a temperature close to the boiling 
point of the fluid and at the same time surface layer 
remains at the ambient temperature. This concept of 
solar ponds provides an exceptionally large area of 
collection of the order of square kilometers at an 
average annual efficiency of 20%. The present collector 
type of solar v;ater heater has an efficiency of 50%. 

The difference lies in the order of magnitude of 
collectinc area and cost. 
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This hind of solar pond is also known as the 
non-convective solar pond, because of the formation 
of a non-convective zone, A solar pond essentially 
consists of three zones as shown in Fig. (1.2). A 
relatively thin zone (UCZ-Upper Convective Zone) 
or surface zone at the top, which contains very little 
amount of salt. The convection is introduced in this 
zone mainly due to wind contact at the surface. The 
second zone, which is the most important part of solar 
pond, is the non-convective zone (KCZ-Non-Convective 
Zone) or gradient zone. As the name implies, there 
is no convection in this zone. MCZ transfers heat 
only by conduction. The third zone is at the bottom 
which is convective and homogeneous (LCZ-Lower 
Convective Zone) . This zone stores energ^^. LCZ is 
present due to the necessity of maintaining NC2. 

Suppose the v/hole pond is non-convective. On a warm 
sunny day, the temperature at the bottom, of the pond 
will be sufficiently high. The salt layers in contact 
with the bottom will be heated to the extent that 
their density is lower than the layers just above them. 
The bottom, layers will therefore move upward. As a 
result convection starts but remains confined to a 
certain thickness only. It does not spread upward 
because above a certain thickness, the negative 
gradient setup by the thermal currents will be off-set 



4. rv. 
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by the possitive gradient of the dissolved salt 
density or concentration. Thus the magnitude of LCZ 
depends upon the amount of insolation reaching the 
bottom. Of course, the formation of LCZ can be 
prevented by increasing the salt gradient, but this 
v:ill be uneconomical as large quantity of salt will be 
required. LCZ is also necessary'- for heat extraction 
purrjoses. UCZ is the principal source of heat loss 
to the surroundings in a solar pond and the convection 
in this zone is mainly due to atmospheric wind. Since 
most of heat is lost from this zone, it is suggested 
to insulate the ponds against this loss by spreading 
transparent plastic covers over the pond surface. 

The zone boundaries, in general, are not 
stationary in time. These move, either decreasing or 
increasing the zone thickness. This motion is usually 
very slow, typically of the order of a few cms per m.onth . 
The boundaries move so as to increase the NCZ thickness 
(or decrease the LCZ thickness) when the temperature 
gradient is lo-w and vice versa i .e ., .decrease in the 
NCZ thickness (or increase in LCZ thickness) when the 
temperature gradient is high . 

This is the solar pond which has been studied 
in some depth in this thesis. 



PARTITIONED SOLAR PONDS (PSP) 


In salt gradient ponds it is proposed to separate 
out LCZ from NCZ by placing a horizontal transparent 
plastic membrane at an appropriate depth . The 
requirement of the salt will be much less because 
maximum salt is reguix'ed in the LCZ which is PSP will 
have pure water. It is hoped that PSP will be more 
stable during energ-y extraction, NCZ will not be 
disturbed during fluid removal and injection. A fast 
rate of energ^^ extraction is possible. Insertion of 
more than one partition is also proposed. Another 
partition is advocated to separate NCZ and UCZ , PSP 
appears to be attractive but larger collection areas 
are not feasible v/ith it. 

VISCOSITY STABILIZED PONDS 


Convection can be suppressed if the viscosity of the 
water can be increased to such an extent that movement 
of water is not piossible. This can be done by adding 
certain gelling and thickening agents to the water. 

The pond contents are now semi -solids. Materials 
suitable for viscosity stabilised ponds should have 
high transmittance for solar radiation, high 
thickening efficiency and should be able to perform 
at temperatures upto 60 °C. It is certain that the 
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proposal will not: be economically competitive with 
salt gradient ponds. Carboxy-vinyl polymer and 
cellulose derivatives have been found useful as 
thickening agents. 

MEMBRANE STRATIFIED SO LAR PO ND (MS SP) 

In flat-plate solar collector, honey comb structure 
is used to suppress convection. Similar principle 
has been proposed in KSSF . The convection in the 
upper layers is suppressed by placing a suitable 
number of vertical miembrane or horizontal sheets or 
vertical tubes [6], The pond now hes only tv70 zones 
LCZ and NCZ . In MSS? no salt is needed anyv-vhere 
making the concept economdcally attractive. 

Maintaining the salt gradient in a SGSP usually 
requires vigilance, and in miany locations salt 
contamdnation is potentially a large environmental 
hazard due to persistence of positive ions in soils. 

In addition to the toxicity to vegetation caused by 
overly salty soils, a major danger of salt pollution 
involves the detrimental effect of salt in drinking 
water upon people susceptible to hypertension and other 
salt related diseases. This concept will be more 
efficient than SGSP due to the absence of UCZ v/hich 
loses a lot of heat to the surroundings , MSSP 
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membrances can be placed very close to the s\irface, 
almost completely eliminating the UCZ to the 
surroundings , Three binds of membranes suggested 
are; horizontal sheet, vertical tube and vertical 
sheet. Teflon is suggested to be the suitable 
membrane material because of its long life, high 
transparency, inert nature virtually to all chemicals, 
commercial availability in all sizes and thicknesses. 
Water is the most preferred liquid but others like 
ethanol, combination of water and ethanol can also 
be used . 

(6) SATURATED SOLAR PONDS 

In this concept, the salt density gradient is 
maintained by the temperature gradient. A salt is 
chosen for which solubility is a function of temperature, 
i The pond water is kept saturated with such a salt at 
all levels and since the pond is hotter at the bottom 
than the top, more salt is dissolved in the bottom 
than the top. As there is no vertical diffusion of 
the salt, S3P requires no maintenance. Kg CI 2 is the 
most preferred salt for this kind of concept ^ Borax 
and KNO^ sre other salts . Such ponds have been built 
on small laboratory scales only. The possible 
precipitation at the bottom of the pond of white salt 
cr^’-stals of MgCl^ will reflect more solar radiation 
rather than absorbing it . 
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1.6 SOLAR POND APPLICATIONS 

Solar ponds collect and store energy for use at 
temperatures belovj the boiling point of the lower zone brine. 
They are lower in cost than the conventional flat plate colle- 
ctors in almost all locations and much lower in favourable 
locations. Typically averaged energy gain per unit area of 
pond is comparable to that of a flat plate collector for the 
same average values of solar input and output temperature and 
the pond advantages are to be found in the lower cost and the 
large intrinsic heat storage capacity. Pond limitations are 
the impossibility of mounting on roofs and the heat losses to 
the earth that make uninsulated small ponds inefficient. On 
account of these losses^ ponds of only a few hundred square 
meters area for individual house or water heating are not 
economically competitive. 

Within these limitations, solar ponds have a wide 
range of potential application for providing low^ temperature 
heat at a low cast . Following are some of the areas where 
solar ponds have been used advantageously. 

(a) Power production 

(b) Space heating 



(c) Crop dry i no 

(d) Domestic water or swimming pool heating 

(e) Desalination 

(f) Industrial process heat 
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CHAPTER 2 


REVIEW OF LITERATURE 


2,1 CONSPECTUS 

UnliKe other thermal collectors, the salinity gradient 
solar pond is a natural phenomenon. It is found in a few 
places in the world vlnere a rare combinerion of circumstances 
results in a body of water more saline at some depth belov7 
the surface than it is at the surface. The first study of 
natural solar heated salinity-gradient lakes v;as that publishec 
by Kalecsinsky [7]. He collected temperat’ure and salinity 
data from several lalces in Transylvania (Hungar^O and 
demonstrated by alalysis and laboratory’- experiments that the 
heating, to temperature as high as 70 ®C, was caused by solar 
energs^ penetration into vrater that Vvas maintained non-convecti\ 
by the density gradient resulting fromi dissolved salt. It is 
interesting that even at that time Fialecsinsky proposed the use 
of artificial solar pond to collect and store solar energy'' for 
domestic and industrial applications. 

Following Kalecsinsky ' s work there was some continuing 
study of the nat’jral solar heated lakes in Transylvania, but 
there was no immediate attempt to act on the proposal to 
construct a solar pond for practical utilisation of solar enerc 
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The SGSP was again proposed as a device for collecting 
and storing solar energy in 1948 by R . Bloch in Israel and 
solar pond research was begun there ten years later by Tabor. 
Bloch conceived the use of SGSP on the basis of the heating 
observed in the Transylvania lakes, first studied by 
Kalecsinsky, This first period of solar pone research in 
Israel continued until 1966 and resulted in very important 
contribution to o'-ir understanding of the problems involved 
in operating solar ponds and in developing their possible 
applications . The primary- goal V7as poxeer production, but 
some studies v;ere made of other possible applications by 
Tabor and Matz [s]. 

Solar ponds were studied subsequently in Australia 
by Davey (1968) [9j, in Chile by Kirschmann (1970) [lO], in 
India by Jain (1973) [H]/ and in the USSR by Usmanov (1971) 

Sc Eliseev (1973) [12], The Australian pond, which was 
operated for two years , constituted a preliminary feasibility 
demonstration. Kirschmann explored the possibility of using 
ponds for electric pov^er generation, process heat, and 
desalination of x\-ater. The work in the USSR consisted of 
theoritical studies and small laboratoiu' experiments . None 
of these early works outside of Israel were continued for more 
than a few years. 

The present period of pond research began in 1974, with 
a theoretical study of ponds for space heating carried out by 



19 


Rabl and Nielsen (1975) [13 ], small pond experiments in Canada 
by Saulnier et al (1975), [14], possible application of solar 
ponds to house es well as food and paper process heating by 
Styris et.al [15], a scheme for modifying the dead sea into a 
huge solar pond by Assaf [16], the performance of the ground 
storage beneath a solar pond by Akbarzadeh and Ahmaoi [17], 

Pond v'ork has been continuous in several sites since that 
time and thei'e is nov; activity in many countries. 

2 .2 ANALYTICAL APPROAOiES 

There have been a number of efforts to solve the heat 
and mass diffusion equations, applied to solar ponds, 
analytically and. predict their performance. A few important 
contributions in this direction are discussed below, 

2.2.1 WEINBERGER'S APPROACH 

Weinberger [iS] v;as the first to study the physics 
of the problem. He developed expressions for the absorption 
of solar radiation, temperature rise, rate of energy v.uthdrawl, 
efficiency and stability conditions of the solar pond . In his 
analysis Weinberger made use of the following assumptions; 

(1) The pond is of infinite size; edge losses are 
negligible and one dimensional conduction holds. 

(2) Thermal properties of salt water and ground are the 


same and constant 
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(3) Natural thermal convection is prevented due to high 
salt gradient. 

(4) No distinction is made between direct and diffuse 
radiation. 

(5) Effect of cloudiness is neglected and 

(6) Salt solution remains transparent even after long 
exposures . 

Based on the above assumptions, the temperature 
distribution T (x,t) in the pond containing internal source 
of heat H (x,t) is governed by the following diffusion 

ecruat-ion 

j|- (K = PC ( - H (x,t) 

Weinberger assumed K, P and C to be constant, hence 

d^T _ 1 dT K (x,t) 

6^2 a " K 

hJhere, 

a = (Thermal diffusivity) 

Weinberger integrated the above equation for the 
following five boundary conditions to predict the temperature 
profile T (x,t) of the pond: 

(1) The solar flux 0 (x,t) that reaches the bottom is 

completely absorbed there . 
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H (x,t) = “ [ ^ (x/t) ] is the solar energy 

absorbed at a depth x in the water body. 

Energy is removed at a rate U(t) from the bottom 
of the pond. 

The temperature of the pond surface is equal to the 
ambient temperature. 

The initial temperature of the pond is kno^-.’n. 

RABL AND NIELS E N'S MODEL 

Ari Rabl and Carl E, Nielsen [l3j made assumptions 
that V7ere very similar to those of Weinberger but were more 
realistic in nature to the problems of solar ponds. They 
considered the presence of a convective zone in the pond and 
also took different thermal properties for the ground and the 
salt water. They thus altered the assumption numbers (2) and 
(3) of the Weinberger approach. They also restricted their 
attention to the steady state condition because solar ponds 
are primarily meant for long term energy,’’ storage. Their 
analysis can be broadly divided into four steps: 

(1) Calculation of time independent component of 

temperature 

For this they solved the steady state heat conduction 
equation and to applied energ-y balance equation to the 
insulating layer (NCE) . 


( 2 ) 

(3) 

(4) 

(5) 

2 . 2.2 
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(2) Calculation of the S ini^soidally varying component 

of temperature : 

This requires the heat balance equation of the 
convective zone to include the heat exchange from the 
ground as well as the heating and cooling of water 
in the convective zone , 

(3 ) Approach to steady state : 

They considered tv/c different cases; (a) slab of 
finite thickness and (b) the semi -inf ini te solid. 

They assumed zero temperature throughout the ground 
and “he surface temperature of the UC layer to be 
constant, and equal to the ambient temperature. 

{ 4 ) Estimation of edge losses in a finite pond ; 

The pond vjas assumed to be circular v;ith radius ' a ' 

and depth 't‘ , The thermal properties of the ground 

and salt vc^ater were ass'omed to be constant. Presence 

of convective zone and time variation were nealected. 

« — 

All radiations were ass-ijmed to be absorbed at the 
bottom of the pond. 

2.2.3 KOOI'S APPROACH 

C.F. Kooi [iP] analysed the therm, al behaviour of a 
solar pond treating it as a steady state flat plate collector. 
According to this model, a solar pond has three zones: 
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(1) UCZ (Upper Convective Zone) having low and uniform 
salt content ana a constant temperature close to the 
ambient temperature . 

(2) UCZ followed by NCZ (Non-Convective Zone) in which 
the salt content increases with depth. 

(3) MCZ followed by LCZ (Lov/er Convective Zone) . 

Here salt content and temperature are nearly constanr 
Kooi, ho'wever, does not consider the edge and the bottom 
losses in his analysis, 

2,3.4 HAWALDER AND BRINKWORTH APPROACH 

Hawalcer and Brinhv.-crth [20l studied the effect of: 

(1) Extinction coefficient of solar radiation 

(2) Thickness of storage zone 

(3) Bottom insulation 

(4) Load application on the performiance of a solar pond. 

They combined the surface heat loss by convection, 
radiation and evaporation into a single term, as: 



n 

'"LS 


- ®a> 


i-fhere 

'^LS 

S\irface 

heat 

loss 



Surface 

heat 

loss coefficient 


e 

s 

Surface 

temperature 


3 . 


Ambient temperature 
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For determining they used empirical equations obtained 

from the measurement of losses from a body of convecting 
fluid at a nearly uniform temperature above the ambient, 

2 .3 FIRST EXPERIMENTAL POND 

The firsc experi'.mrntal poxid vjas constructed by 
Prof Tabor [21] in 1958 in Israel. The results obtained, 
however, v;ere not enccur aging. This v:3s attributed to lov; 
solubility of the corr-mon salt used. The follov’ing aspects 
v.-ere studied in detail: 

1. Amount of solar radiation reaching the bottom 
of the pond . 

2. Condit-ions for stability. 

3. Rate of diffusion of salt. 

4. Extraction of heat frcmi the lower layer. 

5. Effect of v:ind induced waves. 

6 . Effect of evaporation and 

7. Effect of dirt and dust falling on the pond. 

Based on the experiences gained during the project 
the following were suggested as the major problems to be 
tachled, to make the use of solar ponds economical; 

1, Extracting heat from the bottom v/ithout disturbing 

the density gradient . 
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2, Suppression of upper zone mixing due to wind. 

5. Method of keeping the pond clear. 

4. Suppression of biological growth in the pond with 

time . 


2 ,4 PRESENT WORK 


In this v7ork a computer model of a large solar pond to 
be located in the northern pert of India has been developed. 
The input to the program are the hourly solar radiation, 
ambienr temperature and relative humidity data for Mew Delhi , 
Also informations about the depth of pond, initial temperature 
of pond, rherTrial properties of salt water and ground are 
provided . The temperature rise of the storage zone of the 
pond for various rates of energy removal is the output from 
the program. Basic one dimensional heat transfer equations 
are developed incorporatinc the individual terms for ground 


loss, and convection, radiation and evaporation heat loss 
from the pond surface. These are then solved numerically 
■with appropriate initial and boundary?- conditions. An improved 
equation for transmdssicn of solar radiation as suggested by 
P.T. Tsilingris has been used. The effect of the clarity of 
pond’s water is uaken into account with the introduction of 
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a v/ater quality factor. Pond's pcrf crrcnce under various 
rates of heat removal, different load application and 
removal time schedules, variable and constant transmission 
co-efficients and consecutive cloudy days is analysed ana 


discussed 



27 


CHAPTER 3 


THE MATHEMATICAL FORMULATION 


The salt gradient solar pond, considered in the 
present work, consists of three distnct zones: an upper 
convective zone (UCZ-) , a non-convective concentration 
gradient zone (KCZ) , and a lower convective zone (LC2) or 
energy storage zone . The thermal performance of such a pond 
depends on the thickness of each of the above three zones. 

3.1 DEVELOPMENT OF THE MODEL 


The pond is considered to be one dimensional with 
respect to heat transfer, with constant depth ana of infinite 
area. The effect of side boundaries and edges of the pond 
is ignored. This assumption may not cause significant error 
in simulation results because for practical purposes pond's 
lateral dimensions ax'e very large with respect to its depth. 
The model considers the transmission of the incident solar 
radiation through the surface, its attenuation as it passes 
the subsequent w'ater layers and the final absorption in the 
LCZ . The bottom of the pond is assumed to be perfectly 
absorbing and the radiation entering the LCZ is completely 

absorbed there, either in the vjater body or at the bottom. 

Xt is further assumed that the stabixity ox tne pond 

maintained even at high LCZ temperatures. 
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The temperature of the LCZ being higher then that of other 
layers, it loses heat to ground and KCZ. The UCZ loses heat 
to environment by convection, radiation and evaporation. 

In a subsequent paragraph the above mentioned 
processes are discussed in full detail and presented in a 
mathematical formi, using basic principles of heat transfer. 

3 .2 ATTENUATION OF SOLAR RADIATION IN A SOLAR POND 

The performance of a solar pond largely depends on the 
amount of the radiamion vhnich reaches rhe bottom region. The 
UCZ absorbs a large portion of the lone .-.’ave solar radiation. 
This portion of the solar iniaut can not be utilized and is 
eventually dissipated to the atmosphere. 

The most v/idely known mathematical expression on 
transmission of solar radiation through a solar pond in 
literature is the Rabl and Nielsen (R.N.) [l3] four term 
exponential expression given by 

4 

H (x) = X H 2 jj.. . exp (- n. x) (3 .1) 

i=l ^ “ 

VJhere and n^^ are empirical coefficients and have different 
values for different radiation vrave length bands. The above 
expression has been reduced to a simple algebraic equation 

h (x) = 0.36 - 0.08 In x^lO^m x 0 .01 m (3.2) 

by Bryant ana Colbeck [22]. This relation gives the 
per unit attenuarion of solar radiation with depth x of the 
pond. The expression, however, is net very accurate as it 
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has been derived from results based on old and early pond 
data. More recently P.T. Tsilingris [23] has suggested the 
follov.-'ing expression using n'umerical techniques and dividing 
the solar spectrum into 19 non-equally spaced bends 

19 

h (x) = Z exp ( - n,. x) (3,3) 

i=l 

This expression may be fitted by polynomials with a high 
accuracy such as the one of fourth degree, 

h (x) = 0 .67031 - 0.35170 x + 0.19785 x^ 

- 0.05567 x^ 4- 0.0058 x"^ (3.4) 

and this can be fiurther approximated v;ith good accuracy by 
the simple algebraic relationship 

h (x) = a - b In x^ 5.5m > x 1^0.2 m (3.5) 

Vfnere a = 0.46 & b = 0.0953 

The transmission coefficient can be calculated using 
Fresnel relations [2], as: 


t 


0.5 [ 


Sin^ ( i - r) 
2 

Sin^ ( i -f r) 


2 

Tan ( i - r) 1 

o J (-'.o; 

Tan^' ( i 4- r) 


VTnere r and i are computed using the Snell's la'w i.e.. 


Sin i 
Sin r 


1.33 


(3 .7) 


end the incident angle relation, respectively. 

= Sin A Sin 0 4- Cos a Cos 0 Cos w 


Cos i 


(3 .8) 
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v\T^ere, 

A, the declination angle is given by 


A 


23.45 Sin 


(2T1 


284 -f d V 
365 ^ ' 


(3.9) 


d is the day of the year with value of 1 for first January, 
and the hour angle 

= ^ ( 12 - hr) (3.10) 

If we introduce a dimensionless water quality factor, 
Vv^,, [231 which represents for claiitv status or the solar pond 

v;ater, being assumed as 1 for distilled or ideally clear 
natural water, then, eqn . (3.5) can be rewritten as: 


h (x) = 


( a - b In X ) 


( 3 . 11 ) 


Vvhere Vv 


Qf 


= 0.7 to 0.75 for practically pure pond water. 


The amount of radiation reachinc deoth x of the oond 


IS then given oy 


H (x) = h (x) t H 


( 3 . 12 ) 


3 ,3 ENERGY FLUXES WITHIN AND AT THE BOUNDARIES OF THE POND 


Heat energy is lost from a pond primarily by convection, 
evaporation and long wave radiation from the surface of the 
pond to the sky and by conduction to the earth at the bottom. 
Besides, a fraction of the solar radiation falling on the pond 
is reflected back and hence is not available to the pond for 
utilization . 
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There is transfer of heat energy due to interaction 
between the layers of the pond, as they are at different 
temperacures . The LCZ, owing t '. its higher temperature conducts 
heat to NCZ and to the ground. The KCZ, inturn, conducts heat 
to the UCZ which loses it to the ambient air. These energp’- 
fl'uxes are shown in Fig. 3.1. 

In order to simulate the behaviour of the pond, energy 
balances of the various zones in the pond are considered and 
the resulting equations solved simultaneously. 

3.3.1 ENERGY BALANCE OF THE UCZ 

The energy' exchanges of the UCZ with ambient air and 

one NCZ are shcv/n in Fig. 3.2. The UCZ is assumed to be at 

a uniform temperature T . The UCZ receives eneroy* fro.~ the 

u 

sun. Some part of incident energy^ is reflected bacx. The 
fraction of incident energy able to pass the surface layer is 
dependent on the transmission coefficient . 

Hence, 

Incident radiation energy transmitted = t (3.13a) 
Radiation energy^ transmitted through 

UCZ. = Th' h, (3.13b) 

* S J-]_ 

(3.14) 

x=l. 

Energy lost by' convection from the surface 

of the pond, q = h (T, - T ) (3.15) 

C c ^ ^ 


Energy conducted from 


to UCZ 




Non Convective' Zone 


4 Si 


Lower Convective Zone 


Fig 3,1 Energy Fluxes In Solar Poi 


. Depth Dimensions are in meters) 


(All 
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VJhere the wind convective heat transfer coefficient [2], 

is Given bv 


5.7 + 3.8 vh K 


(3 .16) 


Vvhere V is wind velocity in m/s 


Energy lost by radiation, q = £ a F (T ^ x 

i\ s XI s ley / 


(3 .17) 


T - 6 .( 


-g^- lost by evaporation [25l, ^ (3.18) 

“* S V qU X#0U ir’. S a. 

S L. 


fnere [26] f = Exp (18.403 


) , mm Hg 


(3 .19) 


RH .F^ , mrr! Hg 


50ve expressions of eqns (3.13) to (3.19), we can 


write the tonal invvard fi'o:': of UCz, as 


T Ks + K 


ind the total outward flux 


h (T - T ) + e 0 F (T ■ 
c u a s u 


•iky ^ ^1. 


^ 1.6 

S u 


t - r 


The difference in the expressions a & b is responsible 
for the heat accumulation in the UCZ and hence its temperature 
rise. Hence \ve can write 


k + K III T (T 


c u a 


e 0 (T - T , ) 

s u SJCV 


/V h 

H — 5“ (P^-P^)^ (3.20) 
s a- 

s t; 
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by re-arranging and transposing, we get 


c u 
b t 


^ T Hs (1 - hi ) t f 3^ 


a 

kt: 




(T. - a\) - 


a dT 


x=l 


a e a F3 (T^ - T 


K 1, 

X 


'sky 


a (F - F^) 

C Sc 


1.6 C 


K 1 


( 3.21 


\nere a 


K 


, x'nerjTiax cixxrxisxvx'cv 01. sal^ v/arer Cx 


3.3.2 ENERGY BALAI^CS OF NCZ 

The NCZ in a solar pond is a very vital portion, as 
it acts as the insulation and hence reduces the loss of energ-y 
from the LCZ . But as the thickness of NCZ increases, the 
energy flux reaching the LCZ decreases. So for better 
performance of the pond, the thickness of the NCZ must be 
optimised for better insulation and lower attenuation. The 
energy flxix terms in a typical NCZ are as shown in Fig. 3,3. 

The NCZ being at a higher temperature than the UCZ 
conducts heat to UCZ. Similarly because of the higher 
temperature of the LCZ, the NCZ receives energy by conduction 
from; LCZ. NCZ also absorbs a portion of the incomdng solar 
radiation and transmits the rest to the LCZ. For the purpose 
of analy=-sis NCZ is divided into a number of smaller elements. 
The above mentioned energyr flux terms for a typical element 
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in the NC^ can be expressed mathematically as; 


Energy lost to the upper layer = K ^ I ^ 


Energy gained from the bottom layer = K err 

o 


(3.23) 


X-i-AX 


^nerov inout due to sol 


sr raoiatron, using equation 


(3.12), H (x) = h. -crT 


( Z: ^ A \ 


Energy leaving the elements due to transmission of radiation 


‘'I 


(3.25) 


'The difference betVv'een the input and the output terms of the 
elements is responsible for the temperature rise of the element 
Kence, 


PCax-^ = 
he 


dx X-rA>: 


^sJ - [ 


K 


:x 


or 


P c 




^Xi-AX "s ^ 


K 


X I X4-AX 


- K 


dT 


(3*26) 


ZS (3.27) 


or 

PC ^ 
c 


A ^ 

"■'X 


T K ( - ^) 

S OX 


(3 .28) 


But from. Eqn. (3.5) h = a - b In x 


hence , 

dh 

dx 


X 


(3 .29) 
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or 


Substituting eqn (3.29) in eqn. (3.28), we get 




dt 


K 


52^ 

d72 


X 


a.T 

-%t 


a 


2 

d T 
A 2 


xH^b 

_i_ o: s 

K X 


(3.30) 

(3.31) 


The variation of 
obtained by solving eqn. 
conditions applicable to 


ter.perature in the IICZ can then be 
(3.31) subject to the boundary' 
the hCS. 


3.3.3 ENERGY BALANCE OF THE LCZ 


The various energ-y flux terms involved in the LCZ 
are as shown in Fig. 3.4. Due to presence of convection, this 
zone is assumed to have a constant temperature T^ . In 
mathematical form these interactions of energy fluxes can be 
written as follovk's; 


Energy?" conducted to the NCZ, 


Energy conducted to the ground. 


= K 


5T 

A X 


x=l. 


•q 


= K. 


6 X 


x=l. 


3 

Energy extracted for useful applications = 
Energy received by the LCZ = h-, ^ 

Hence, total inward flux = h, T H’ 

io s 


(3.32) 


(3.3 


•^3) 


( o .d 4) 

(3 .35) 
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T 

and total outv;ard flux = K ^ _ . + K_ - - ? • — I , + q 


(3.36) 


The difference between expressions (3.35) and (3.36) is 
responsible for the heat accumulation in the L-CZ, hence, 

(I3 - Ij) PC - [k|| q^: 


(3.37) 


3y transposing v;e get. 


d^L 
g t 


h, T H 

^2 K . ^ j 

PC (1. - ± 2 ) ~ (1. - l.)PC ‘ arx j x=l., 




PC (1„ - 1.) 


o G 
dx 


PC (1^-1^) 

O ^ 


(3 .38) 


Sqn. (3.38) can then be solved in conjunction v/ith 
other equations subject to appropriate boundary conditions to 
predict the behaviour of the LCZ. 


3.3.4 ENERGY BALANCE OF THE GROUND LAYERS 


The ground layers receives energy only from the LCZ 
and conduct it to their lov/er layers . Since there are 
temperature variations with depth as well as time, the ground 
is divided into smaller elements parallel to the lo'wer 
convective zone where in temperature can be assumed to be 
uniform with depth without incurring appreciable error. The 
energy flux terms for such a layer are shown in Fig. 3.5 and 
can be vj-ritten mathematically as follows: 
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Inerg^r conducted from the upper ground layer 

rri 

I 


K. 


dx lx 


(3 .39) 


Energy conducted to the lower ground layer 


- I X+AX 

The difference between (3.39) and (3.40) is then 
ecuated to the rate of enercr/ storace 


(3 .“fiO ) 


pc 

- 

CD 

1 





C X 

or 





pc 

d^G 

^ Tt = 


r 

L 

0 

d X 


X 


(-K, 




aT 


X+AX 


3y transposing we get. 


G 

33c 


X-rAX 


X 


d'^G 
d t 


_2 G 


6x^ 




(3.41) 


(3.42) 


(3.43) 


i'ihere = 7 ^ • Thermal diffusivity of the ground 

Cjr 

Eqn. (3.43) simulates the temperature variations of the ground 
lavers . 


3.4 


BOUNDARY CONDITIONS 


Temperature distribution in the solar pond is then 
obtained by solving the following equations simultaneously 
subject to the boundart’ conditions listed below. 
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o u 
d t 


^ TH^ ( 1 - h, ) + f 


1, OX 
1 


X = 1 


kl. 


(T - T ) 
u a 


a 


a (T 


u 


kli 


) 

sky 


a k (ir'g - ir-_) 

___ _______ 

i Sr 


(3.21) 


dT 

t 




£ 
■ K 


T H^b 


3 .31) 


>T- 

$ jj 

d t 


u T H 

i-. S 

z 


PC (1, - 1.) 


K a T 

(l^ -l„)pc 'x 


x=l. 


ana 


513 

at 


C^c 


pc ( 1 . - 1-) :■ X I x=i„ pc {!-- 1^) 


0.r 


i2 T 
6 -^G 


(3 .3S) 


(3.43) 


( 1 ) 


The upper convective zone surface temperature is 
equal to ambient temperature at ell times i»e., 

rp _ rr. _ S p 


1 ,44) 


(2) The temperature of the KC2 at x = 1^ is equal to the 

temperature of the UC2 at ell times i.e,, there is nc 
abrupt change of temperature at the interface. 

This can be written mathematically as: 


T (x) = T^ at X = 1^, t ^ 0 


(3.4 5) 
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(3) Temperatxire of the KCZ at x = 1-, is the same as 
the temperature of the LCZ 

T (x) = T. at X = , t > 0 (3.46) 

Li <L 

(4) The bottom of the pond is assumed to be at a 

t:6rnp0rat:‘ure er*ual "co that r.f the LCZ, 


( 5 


/ 


(x) 


JL _ O tU- 




\-J ♦ tr / / 


Temperature of the ground belovj a certain depth from 
the bottom of the pond (1.0 m in this case) is equal 
to the ambient temoerature . 



1- + 1.0/ L 0 


( 6 ) 


The pond is assumed to be at the ambient temiperature 
throughout its depth initially/ i.e.. 


u 

T (x) 


T^ (x) = 

<J 


¥ X 


¥ X 


Y 

n 

3t 

JC 


t = 0 


(3.49) 


Thus for the solution of the problem eqns . (3.21)/ 
(3.31)/ (3.38) and (3.43) have to be solved 


simultaneously subject to the boundary’' conditions 
given by eqns. (3.44), (3.45)/ (3.36), (3.37)/ (3 .38) 

and (3 .49) . 
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3.5 NUMERICAL MODELLING 

Follov;ing relations are used to non-dimensionalise 
the governing equations. 


X® 


X 

^3 

Non-dimensional 

depth 

1 


= 

at 

.2 

^3 

Non-dimensional 

time 



=: 

KT 

H i., 
a 3 

Non-dimensional 

temperature 

(3.50) 










Non-dimensional 

solar 




a 

radiation 



po 

= 

P 

^t 

Non-dimensional 

pre s sure 



Superscript ‘ ® * is used to denote the ncn-dimensional 
forin. On substituting the above relations of eqn. (3.5C) 
into eqns . (3.21), (3.31), (3.38) and (3.43), and simplifying, 

the respective governing equations take the following fonr;. 


aT « 
0 u 


1 at= 


c T° 


h 1 

tK^" (1 -h^^) (V -T^o) 


^ > o 4. ^ A 

- P (T - T . ■) - Y (P ° - P ") 

u SKy s a 


(3 .51) 


6t° 


dt° 

(1 - 2,) 


~L° 
2^ d t® 


n 2 H ® 

d T® T - s 

^ n2 X® 

Cx° 


V, _ tj o _ " 

^2 ^ “s 6x° 




^''G d-^G® 


1.0 


(3.52) 


- (3.53) 
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ttp 




( 3 . 54 ) 


<here , 


I1/I3 


( 3 . 55 a) 


I-7/I' 


e ^ s ^3 


(3 . 55 c) 


Similarly on non-dimensionalising, th 


a oounuarv 


:onaj.tions can oe vvrzLtten as: 


for ail 11 ° ^ 0 


/v^ — 

- 


x° = 2 . , t° X 0 


T° Cx) 


=t x° = Z^, t° X 


it x° = 1 . 0 , t° >/ 


^ , t° ^ 0 




,.6 C H, 


( 5 . 55 e) 
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at t° = 0 


(3 .61) 


It is ver-'/ tediou 




xme consuminc 


so xve 


governing equations analytically. However, their finite 
difference analogues can be solved easily. For v.'riting 1 
finite difference equations Crank - Nicolsen [27,28,29] 
scheme has been utilised. The governing equations qhen 


assume tne rorm as erven belovj. The 


■script ' ° ' usea 


to represenp the non-dimensional form is being droped for 
convenience . 


t (#^) r ^ 


+ t'' (1 - h, ) 

s 1^. 


-£ £ (-T ^ _ .y _ p (g _ g 

K u a u skv 


- T (P3" - ] 


(3 .62) 


n+1 

-r T, , + 2 (1 e r) t/ 

X “"-L 1 
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n-f 1 


'L 


T ^ — ) r h 

•^L ^1 - L 


1. ^ 


n _ n 
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CT 
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*T 

O-^G 


K "5^ 


n 


- T) 


1.0 


(3.64) 


•o 1 ^ 


-1 


(1 


s) T. 


n-f 1 n-f 1 

- s T i . , = s T 


n 


j-1 


2 (1 - s) T." + s 1.^;^ 


(3 .65) 


iere r 


^t/Ax' 


V. ♦ O O o. / 


a.. 


:^:) ' 


/ o, 

\ o « c ■>«' xv J 


Superscript 'n' depicts the time level and subscript 
'i' end 'J' represent the space level in the NC2 and in the 
ground, respectively^. The above finite different equations 
represent the numerical model of the solar pond. 


3 .6 THE COMPUTER PR0GRAI4 

Eqns (3 .62), (3.63), (3 .64) and (3.65) vvere numerically 

solved on the ccmputeo, HP 900Q/S50 system. The follov;ing 
paragraphs are devoted to the explanation of the development 
of the program. 

The model employs a thermal netvvork that divides 
the solar pond into horioontal sublayers. The KC2 is divided 
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into sublayers each v.dth de^jch of 0.1 mv;hile the UCZ and LC2- 
are one sublayer each as there is no variation of temperature 
with depth in the UCZ and the LCZ. Network nodes are located 
in the middle of each sublayer and the node for the ambient 
temperature is at the top of the pond. Nodes of adjactent 
sublayers are in conductive triernel contact and -he temperature 
variation betv>reen nodes is assu-med to be linear. Physical 
parameters are assirmed to be constant throughout each sublayer. 

Since hourly teirperature , solar insolation and FC-I 
data for New Delhi were available, the time interval for 
iteration was taken as one hour. The computer was fed with 
one full year solar insolation dara beginning from January^ 1st, 

The pond depth, thickness of the NCZ, the thickness 
of the UCZ, and thermophysical properties of the salt solution 
and the ground, were given as input to the program. 

The first fev; steps of the prograrri have been devoted 
uc the conversion of the variables into non-dimensional form 
and development of tria-diagonal matrices required for the 
Eqn. (3.63) and (3,65). Calculation of the ancle of incidence 
using the latitude of the location, time of the day and day of 
the year plus the computation of the coefficient of transrrdssion 
using fresnel's relations have also been included in the program. 
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The tria-diagonal matrices are decomposed into the upper and 
the lower triangular miatrices. Equns . (3 .62), (3.63), (3,64) 

and (3.65) are then solved with iterated value of gradients. 
The output of the program is obtained in the form of the 
temperature of LCZ at every 350 hour internal. The program 
v^as run for a number of heat exuracuion loads on the pond 
with variable davs for the application and removal of loads. 
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CHAPTER 4 


RESULTS AND DISCUSSIONS 


4.1 EFFECT OF TRAI^SHISSION COEFFICIENT 

The TTiOdel deve loped in this thesis ccTupares nhe 
effect of transiTiission coefficient on the temperature rise 
of the LCZ by calculating the angle of incidence, hence 
transmission coef f ici3n~-, rt each hour of the year, \^-ith a 
consranr transmission coefficient, corresponding to 2 PM 
solartime at equinox [29 j ( t = 0.979), throughout the year, 
vith a constant transmission coefficient of 0.979, the LCZ 
temperature is slightly on the higher side than with a ^ 

variable transmission coefficient. The ma>:iiriUin temperature 
difference of about 8°C occurs at the end of the year. 

Figs. 4.1(a) and 4.1(b) show the variation of 
transmission coefficient with time. It can be seen from the 
figures that between 10 to 14 hrs of the day, the transmission 
coefficient has nearly a constant value for each month of the 
year varying slichtly over the months. Fig. 4.1(b) also 
displays that t increases upto the month of July each year and 
then starts decreasing. The maximum variation in t , at solar 
noon is from 0.961 to 0.979 in a year. Thus employing a 
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Fla, 4, 1(a): Varifition of for the first vSix months of the vear 
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constant value of T = 0.979 leads to higher values of LCZ 
temperature. This can be seen clearly in Fig. 4.2 vdiic'n 
gives the LCZ temperature variation v'irh time using variable 
and constant transmission coefficients. This is in full 
agreement v.’ith the findings of J.H. Hull [29] that hourly 
variation of transmission coefficient over the year has very? 
little effect on the pond's overall performance. 

4 .2 EFFECT OF EVAPORATION LOSS 

Fig. 4.3 shov;s the effect of evaporation loss from; 
the pond surface. Lov:er cur\?e depicts the LCZ temperature 
variation v.'hen all the losses from the pond i.e., convection, 
radiation, evaporation and ground losses have been considered 
Upper curve gives the LCZ temper at ’ure vjhen all but evaporetio 
losses have been taken into account. It can be seen from; the 
graph that around 4200 hr i.e., during June, evaporation loss 
reaches a maximumi. This explains the dip i.n the lower curee 
despite the fact that solar radiation during this period is 
maxim-um! . As the evaporation loss reduces with timse the 
temperature of the LCZ rises till approximately 7000 hr i.e., 
up to October end, after which it goes dov.'n again as radiatio 
losses incres-se fromi the pond surface due to decreasing 
amibient temipereture coupled v;ith decreasing solar insolation. 
At the end of the year, the temperature difference in the LCZ 
under the above two circumstances is only 10 °C, which is not 
very substantial. Evaporation loss, however, is predomdnant 

mainly for the three months beginning April. 



HEAT OUTPUT = 0.0 

-# — • — normal year with evaporation loss 

200.00 n ^ — normal year without evaporatbn loss 


54 



(9pDj6nu93 ui)'dp)3i 3NOZ 3AI103AN00 HBMOI 


ivf.t.ect crP fiv.iporntj on loss on LCZ temperature. 
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4 .3 EFFECT OF ONE MONTH OF CONSECUTIVE CLOUDY DAYS 

In the present study, response of the solar pond to 
several consecutiva cloudy days has also been investigated. 

This is quite important for the application of solar ponds 
for house heating, process heating and electrical pov^er 
•generstion. It is assuried that July in the Indian suc- 
continent is fully cloudy and s'iar insolation fcr these 31 
days is zerc . The temperaoure rise of the pond is plotted 
fcr such an exceptional year and is compared v.’ith a normal 
year in Fig. 4.4 for the case of no heat extraction from the 
pond . It is observed that at the end of thirty one consecutive 
cloudy days, the temperature of the LC2- reduces to a rrdnimuiTi 
of about 105 °C. This is follov;ed by a period of warming up 
as the consecutive cloudy period ends and solar radiation 
starts coming in. The tem.perature of the LC2 further 
decreases as the winter sets in. At the end of the year, the 
temperature difference of the pond for 'a normal and an 
exceptional year of 31 consecutive cloudy days is less than 5°C, 
From the results of Fig. 4,4 it is concluded that solar ponds 
are relatively reliable solar collector-cum-storage systems, 
since even in an extreme case of thirty one consecutive cloudy 
days, with zero solar energg- influ>r, only slightly reduced 
LCZ temperatures are available at the end of the year. It is 
of course clear that probability of having thirty one 
consecutive cloudy days is almost zero. Even for a cloudy cay 
the solar radiation available is almost 15 to 20 percent of a 


correspondj ng sunny day 




one month of conr.ecutive clondv davn 
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A .4 


EFFECT OF WATER QUALITY OF THE SOLAR POND 


In Fig. 4.5, the effect of pond water clarity on LCZ 
temperature is shown. The upper curve is for the condition 
when for the pond water is unity i.e., ideally clear 

w’ater and the lovrer cunt's for the condition when W , = 0.7 5 

pr f 

i.e., a condition of prectically clear v.^ater at which pond 
could be maintained. Due to dirt, floating organic matter, 
algae growth etc. uhe direct solar radiation component 
penetration into the pond is reduced and the diffused radietic 
component is augrrientec due to partial scatuering of the beami 
component of radiation. At 6500 hrs for uhe w'ater quality 
curve of 0.75, the peak temperature is reduced approximately 
20°C. It is clear from^ the figure that for practical ponds, 
the LCZ temper atirre is appreciably lov;er throughout the year 
than for the ideal water ponds. This suggest that for higher 
LCZ temperatures, a higher operating cost is expected for 
maintaining the clarity of the pond v;ater near about unity. 


4 .5 EFFECT OF HEAT EXTRACTION ON SOLAR POND 

A detailed analysis of solar pond performance under 
variable conditions of load, load application day and load 
remioval day has been carried out. For all the cases, all 
the losses and ideal w^ater quality have been taken into account. 
After an initial warm up period, a fixed percentage of the 
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annual mean solar radiation is assumed to be removed 
continuously from the pond at its bottom. It is assumed that 
heat extraction process dees not produce any instability in 
the pond. 



40, 6 0 , and 80 W / 


removal increase 
Therefore, for a 
removal should b 
maintained at a 
LCZ temperature 
and 60 '■i/m'’' , the 
increase being 1 
by the fact that 
total heat extra 
the 80 Vv/m^ load 
the solar inout 


, depicts the case v.'hen load is applied after 
no up period of sixty days, at the rare of 20, 
m" throughout the year. As the rate of neat 
s, the pone temiperature drops considerably. 

practical pond operation, the rate of heat 
e such that the temperature of the pond is 
reasonable level. At the end of 60 days, the 
reached nearjy 90 °C. ?or the loads of 20, 40, 
LCZ temperature goes on increasing, the 
ess for the higher loads. This is explained 
solar inyut at these loads, is more than the 
ction (inclusive of all losses) . However, for 
, the LCZ temperature does not rise further as 
at this load is less than the total heat loss. 


Fig. 4.7 protrays the case v.^hen LCZ has already 
attained its maximim possible temperature after initial ISO 
days of heating. Mext the solar input starts decreasing v^hile 
the total load is increasing due to lov.-er surrounding 
temperatures. As negative temperature gradient has already 
set in, the increasing load rate enhances the temperature drop 
rate throughout the year. 
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Fig. 4.8 depicts the case v/hen the load application 
day is xhe 180th and load is removed on 300th day of the year, 
so that during the v:inter period pond can regain some last 
temperature . Curves shov/ a rising trend after 7200 hrs because 
due to lovjer temperature of the pond, the losses are red’aced . 
Hence solar irr'ut is more than the losses from; the pond. 

Fig. 4,8 also suggests that ue can expect a higher rate of neat 
removal for a small period still maintaining a practical pond 
temperature . 


Figs. 4.9 and 4.10 show the case of the load app-lication 
day as 240th and the load removal day as 365th and 30Cth day, 
resjpectively . These figures suggest that it is not advantageous 
to increase the day of application of load as the useful days 
of pond's operation (days of heat extraction) narrov; dcrwn with 
out any comipensating advantage. Further a high pond tem.perature 
easily destablises the pond due to starting of boiling v/hicn is 
not permissible in an actual solar pond. It is clear that for 
practical operation of a solar pond, excessive heat must be 
removed as useful heat for various applications. VJater should 
not be allowed to boil since it will destroy the salt gradient 
ana hence the HCZ . 


4,6 


OINCLUSIONS 


Based on the results obtained from; the present 
computer simulation of the performance of solar ponds in the 
northern part of India it appears that solar ponds are reliable 
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Pig. 4.6: LCZ temperature v/ith the loacJ applied on the 60th day and 
load removed on the 365th day . 
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Eio. 4 . 8 : LCZ temperatiirr^ wH:h the load applied on iRoth dev 
removed on Booth dav. 



LOAD APPLIED ON 240th DAY REMOVED ON 365th DAY^^— ^Meat output = 20 

— ^Heat output = 40 

■*~~*~-Heat output = 80 

160.00 n ''^-^Heat output = 100 
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LCZ temperottire with t.ho load applied on 24nth day, 
removed on 365th dav. 



LOAD APPLIED ON 240th DAY REMOVED ON 300th DAY-'^-^Heat output = 120 w/m*m 

^Heat output = 140 w/m*!n 
”* •—Heat output == 1 60 w/rri*rn 
60.00 “1 •—Heat output = 180 w/m*m 
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Pig. 4.10: LCZ temperature with the load applied on ?.40th day, 
removed on 300th dav . 
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solar collectors for production of hot v^ater all year round i 
this area. As can be seen from Fici. 4,2, the temperature at 
the bottom of a 1.5 m deep pond whose operation starts on 
Jan 1, rises rapidly and reaches a value of 150®C in May with 
no load conditions . 


Meat should be extracted from the pond to the extent 
that boiling does no't occur vfoich 'wo’.^ld disturb the salt 
concentration Gradient . 



Variable has little affect on the pond ove: all 
perf onriance in comparison to the constant i . 

The evaporsition loss from; the surface of a solar pond 
is quite significant during the mo nth s of May, June and July. 

Studv of several consecutive cloudy days shews that 
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4 .7 SUGGESTION FOR FURTHER WORK 

The present study is bssically a theoretical 
investigation into the performance of salt gradient solar 
ponds , The next logical step tov.-crds extension of the present 
v.’orh v;:'u1g be tc construct ana study an experimentei pond to 
verify vne results arc also to relax sonia of tne ass urriOtions 
underlying the pu'esent investigation such as consoenu thermal 
conductivity , one dimensional conduction and constant v;ater 
clarity iN'ith time. 
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APPENDIX A 


• thermal stability of solar ponds 

Vveinberger [18] has developed a one dimensional thermal 
r.odel and explained the stability of solar ponds. SP tvill be 
stable, i.e. all the three zones remain present, when the 
density gradient on account of the salt concentration is 
equal or greater than the negative density gradient produced 
by the temperature or the total derivative of density v;iuh 
respect to depth is greater than on equal to 0 . 

P = f (C, T) 


d>c 




^ u. ( 62 ) ^ > 0 

c'T dx ‘ oT c <i>: 


I'ihere , 


QC 

dx 


n 

Integrating, 


pT dT 

B d>-; 

■ s 


(4^/p) , Thermal expansion coefficient 


of waver 


/ dP\ 1 1 4_ ■ ^ . 

"Sc P ' ‘=’=1'*^ expansion coemcienv 


(aC) „ 
min 


AT 


This expression relates the rrdniimim concentration 
difference required to maintain thermal stability betv;een 
tv;o layers of liquid with aT temperature difference. 
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This equation ensures only thermal stability. In 
long run other kind of in stability also occurs. This occurs 
because thermal diffusivity is much greater than molecular 
diffusivity. So that concentration fluctuations will last 
almost 100 times longer than temperature fluctuation. 

Some natural mechanisms also induce convection cells v.Tiich 
grov; with time. Unless all grcv;th is prevenred, ir v;ill 
eventually destroy the r.'CS . To suppress oscillatory 
disturbances the salt concentration gradient snoold fulfil 
tr:6 follo'wing condition: 



'nnere 


'randtl nurnbej 


Gl 

_ __s _ Coefficient of salt diffusivity 

"^s aT Coefficient of thermal diffusivity 

The minimum salt concentration required to prevent 

oscillatory disturbance is therefore 
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(a) Concentrated salt solution should be periodically 
added at the bottom to compensate for the loss of 
salt by diffusion. 

(b) The top of the pond should be frequently flushed to 
maintain a very lov: concentration of salt in top 
layers . This helps in maintaininc the transparency 
of the pond. 
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APPENDIX B 
NOMENCLATURE 


ABSI 

AK 

m:g 

AL 

ALl 

AL2 

ALPHAG 


DSLTITD 

dsl:< 

DEU'-aiD 

FS 

GRADE 1 
GRADE 2 


Ernissivity of the Pond surface 

Conductivity of salt solution 

Conducti\’'i ty of Ground 

T'otel ooptr. of the pond 

Depth of UCZ-NCE interface 

Depth of NCZ-LCZ interface 

Thermal Diffusivity of salt solution 

Thermal Diffusivity of Ground 

Time intem’al 

Non-dimensional time inter'val 
Element t'nichness 
Non-dimensional element thickness 
Shape factor 

Temperature gradient at depth 1^ 
Temperature gradient at depth I 2 


grad: 


o 


KA 

HC 

KS 

I 


Temperature gradient at the bottom 

Annual average radiation 

Convective heat transfer coefficient 

Solar radiation 

Number cf sublayers in NCE 

Total number of days of cpereticn 



Number of sublayers in the ground 


KA 

mA 

PKI 

QE 



Time of operation in Hrs. 
Day of application of load 
Day of removal of load 
Latitude of the nlace 


Hee 

t extrac 

:tion 

rate 




Amb 

ient text 

pe nature 




Non 

-dimensi 

onal 

arffoi 

ent 

: terrpe 

ratur- 

Non 

-dimensi 

onal 

grou 

.nd 

temper 

oture 

Non 

-dimensi 

onal 

T r'9' 

terr 

iperatu 

re 

Non 

-dimensi 

onal 

sky 

terr 

■peratu 

re 

Non 

-dimensi 

onal 

T'T'O, 

ten 

iperatu 

'>^C: 

Non 

-dimensi 

onal 

NCS 

temperetu 

re 



X 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

ccc*: 


Guxae 

Dr. P.N.Kaul 
Ast T . Professor 


c 

CrrCCT or various PARAnCTCftS on the rCRFORnANCC OF c 

SALINITY GREADIENT SOLAR FOND C 

C 

c 

By C 

Yogesh K Gupta C 

H.Tech. (Thermal) C 

C 

c 
c 
c 
c 
c 

rCCCCCTCCCOCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
dimension hs f A QO , 55 ) , ns 1 c Beo 0 i , ta { 4 oo , 55 J , ta I not BBOO > , tsKynd(BBOO J 
dimension txnd{40)>tgnd(40),a1C40,40),bl (40,40),bttC40,5) 
dimension sip(40),c1C40),a11(40,40),a1y(40,40),y1C40), rhsi C 40 ) 
dimension b£(40,40), b£x (40,5) ,c£(40),a2iC40,40), a2u (40,40) 
dimension dsc( BSOO ) , rhs2( 40 ),aEC40,40),y2(40),rh(8©00) 
dimension rh 1(365,24) 
openfaniT*1 ,file='ny6.dat' ) 
open(unit=2>file='sp. out ' ) 
read ( 1 »)al,alt,al2 


nCCMANICAL CNGINCCRING DCFARTNCNT 
INDIAN INSTITUTE OF TECHNOLOGY 
KANFUF - BOB 015 


reed( 1 ,*)hc,ak,akg,alpha,alphag,he 
reed C 1 , * )absi , f s , phi 
reed C 1 , ) i 


read ( 1 . ♦ ) j 

read { 1 , )na, de i x , del t , ibp 
read (1 , ♦= ) nroe , nma 
do 10 11 = 1 , ibp 
do 10 .1 j * 1,24 

hE( i i , 3 3 )= 0.0 

1 0 cont i nue 

do 20 11=1, ibp 

read <1 , »* ) ( hs ( i i , j j ) , j j *6 , 1 9 ) 
20 continue 

do 30 2 i* 1 , ibp 
do 30 j j«l ,24 

jne as£45> ( i i *-1 ) + j j 

hsi ( me ) ®hs ( i 1 , j j ) 
dscdrsa) =h£l(ma)/ha 
30 continue 


do 40 k= 1 , 1 bp 
readC 1 ,*}(ta<k, l),l=t ,243 
40 continue 


do 50 11=1 , ibp 
do 50 j j = 1 , 24 
m=24* ( 1 i- 1 ) + j j 

ta Ind ( m ) = ( t a C i i , j j )*»'273 ) ♦ak/( al ♦ha ) 
tskyndsim) = (ta<ii,jj}+867)^ak/(al»ha) 
so continue 

del xnd=de 1 x/al 

del tnd = de It ♦a Ipha/ ( a 1 ♦a 1 ) 

z1 =ai 1 /al 

z £ = a 1 2 / a 1 

t undo = talnd( 1, > 

1 1 ndc»t a 1 nd ( 1 ) 
dc 60 n=i , G 
T xnd? n)=t alnd( 1 ) 

60 continue 

do 70 n=1 , 1 0 
t gnd( n )=talnd( 1 ) 
com 1 nue 
do 80 ii 1=1, 365 
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ri?adC 1 , ) (rhi { 1 1 i , jj j ), jj j«l ,24) 

60 continue 

re&dC T , )qe 
gradz 1=00 
gradzS=0 , 0 
gradz3= 0 . 0 
s i gma«S . 67e-06 

bet a«abs i *f s^signa’*^ ( al/ak>**4*ha>*'’f3 
r*deltnd/ ( delxnd^delxnd) 

**r^alphag/«lpha 

c h1 = .75C , 46“0. 0953»alog(al1 )) 

ht = ( . 46-0 , 0953»alog<al 1 ) ) 
ha= (0.46-0. 0953*alog( al3) ) 
c h£* . 75( 0 46-0. 0953’»alog(al£) ) 

do SO iT-.= 1 , 8 
a1(jr-.iT»5«2.0»r + S,0 
bt im ,m)^E . 0-S . 0*r 
SO continue 

do too m ®£ , & 
a 1 ( HI)- t , m > * -r 
a 1 ( jr< , «- 1 ) *-r 
b I ( m - 1 , n ) = r 
b 1 C n , m - 1 ) = r 
I 0 0 cont 1 nue 

do I I 0 m=l , 1 0 
a£Cm,m>*£ . 0+2 . o*s 
b2( m , m J =2 . 0-E . 0*s 
T 1 0 cont 1 nue 

do 120 iti«2 ^ j 
a£( m- 1 . fl> ) =-s 
ae t m j, n- 1 ) »=-s 
02 c m - 1 , in ) »s 
be ( m , m- 1 J *=s 
120 continue 

do 111 nn«*2, 2760 
1 nn=nn/e4 
8n=nn 
ai nn= i nn 
f nnssan/S-*! . 0 
rrnn=rnn-einr: 
me= I ■•■inn 
dme=ffie 

do 1 = 24 o^rrnn 
dol 1 =do2 - 1 
nffi=nn- 1 

i day* int ( nn/£4 ? 

1 hr=rjn-£4 ’i' 1 day+ 1 
rh( nm ) = rht ( i day , ihr ).-^t 00 
i f < dol le . 6 . 0 )goTo 618 
i f ( do 2 . ge . 1 8 )got o 618 
pi=3 . 141568 

omega=Cpi/1£.0)>^=(dol-t2.0) 

de 1 1 a = c3 45* e i n < £ O^pi * ( £84 . 0-i-ame )/365 . 0 ) 

d€ 1 1 ar=de It a^pi /1 80 . 0 

phir=p2»«'phi./'180.0 

cos i rj=s i n ( de 1 1 ar ) ♦E i n( phir ) ■♦-cos ( del tar ) *005 ( phi r) *cos( omega ) 

einci=acos(coEin) 

sirjln*siriCainci ) 

sinre=s2nin/l , 33 

rer = asinf. sinre ) 

a5r*a i nci -ref 

apr=a 1 nc i +reT 

Ft ou= 0 . 5** ( £ in ( asF )/s in( apr >)=♦»♦=£*( I . O-*- < cos ( apr )/cos ( asr ) )^»£ ) 

t OU* 1 -FT OU 
C tOU*r.7T7 

go TO Cl 7 



618 t ou»<) .0 

619 n*«nn-t 

i f C «# . It . n»a )go to 130 
if {«« . gt . n»«)go to 130 
•t 4i*q*/ha 
go to 140 
130 •t«*0.0 

140 do eaa id-i,a 

lond««£300 
c»»1 .0£16 
id«y 1 *id«y+l 
ihr I “ihr^S 
tt*«t*( id«y 1 , ihr 1 ) 

pa*rh{n»)*#xp( 18. 403-3885/ (tta-^-SSO) )/760 . 0 
ps*p«/rh(n®) 

ga»4i*l<^»d«*hc/Ch«*CE*760 ) 

tundn*tundo+(deltnd/ 2 l ) * < gr«dz 1 - ( hc**l/jik ) *( tundo-t«lnd(n») )- 
1b«t«»( tundo**4-t5kynd(nBr)**4)+tou*d«c(n»)*( 1 . 0-ht )-g«ii4i*(p»-p«) > 
t lndn»t Indo+C deltnd/( 1 . O-zS) )* ( h£*tou*d«c (n»>-gr«dz£-Cakg/«k ) 

1 ♦gr«dz3-eta ) 
do 150 i j*t , 8 
j i*1 

b1t(ij,ji)*0.0 
do 150 k*1 , 8 

bltCij^jil^bltCij^jil+blCij^kl^txndCk) 

150 continue 

*i*0 . 18*tou*deltnd»dscCn«)/delxnd 
do 160 k«1,8 
a k * k 

c1(k)*0. 0 
sip( k )*si/ak 
160 continue 

cl <1 )»r*(tundo+tundn) 
cl { i ) *r* ( t Indo+t Indn ) 
do 170 i n*= 1 , 8 
ji*1 

rhst (in)®b1t(in,ji)tsip(in5i'Cl(in> 

170 continue 

a1 l( 1 , 1 >*a1 (1,1) 
a1u( 1 , n*1 .0 
do 180 1*£,8 
aU(l,l-1)»«1(l,l-1) 
e1u{l,l)*1 .0 

a1u( 1-1 , l)«a1( 1-1 ,l)/a1 1(1-1 , i-1 ) 

a1 1 ( 1, 1 ) «a1 ( 1, 1 )-«1 1 ( 1, 1-1 )ea1u( 1-1,1) 

180 continue 

y1(1)«rhs1(1)/atl(1,n 
do 190 »*2,8 

y1 (»)«(rhs1 C»>-a1HiB,ii!-1 )»y1(tt-1 ))/a11(«,*) 

190 continue 

t xnd< i )«y 1 ( i ) 

■in*8-1 

do £00 k*1,ain 

»in=8-k 

»p««in+t 

txnd(»in)*y1(»in>-txnd(»p)^«1u<»in,»p> 

£00 continue 

do £10 k* t ,10 
b2t(k,1 )*0.0 
do £10 1*1,10 

b£t ( k , 1 )«b2t ( k , 1 )+b2( k , 1 )*tgnd( 1 ) 

£10 continue 

do 220 k«1 ,10 
c2(k)*0.0 
220 continue 

c£( 1 )*»*(t indo+t Indn) 



3Q 


cS( j )*s>»=Ctaind{nfn)-»'TalndCr»n) ) 
do £30 k=1 , 1 0 
k j = 1 

rhsS ( K )*b£t ( k , k j ) +c2 ( k ) 

£30 continue 

821 ( 1 , ! 3=a£( I , 1 J 

a2u{ 1 , 1 )*1 .0 

do k*2,l0 

851 C k , k~7 >«a2C k , k-l ) 

«2u( k , k 3 = 1 .0 

*eu(k-t ,k>«82Ck-i , k )/a2I C k-1 , k-1 ) 

8£3 ( k , k )«82( k,k)-a21(k,k“l )»«2uC k-1 , k 3 
240 continue 

y£C 1 )*rhs2( 1 )/a2i (1,1) 
do £50 k*£, 10 

v2( k 3 * f msec k )“«21 ( k , k“l 3 **=y2{ k-l 3 3/a21 ( k , k 3 
250 continue 

tgndt j 5 = v£t j 5 
mo= t 0- 1 
do 260 k= 1 , mo 
niffi=l 0“ k 
np»nitri+ 1 

tgndCn 2 m)«y£<nim)--Tgnd(np)» 8 £u(nim,np> 

260 continue 

gr 8 d 2 l=(grad 2 l+(talnd(nm)“txnd<l ) ) /del xnd ) /£ . 0 
gr«d 2 £=( fgrad 2 2+(txnd( i )’-txnd(i-“1 ) )/delxnd>/£.0) 
grad23=(gradz3+(tgnd(£)“tgnd( 1 ) )/delxnd)/£ . 0 
222 continue 

t undo«t undn 
t indo*t men 
i8n=nn/350 . o 
8nn*nn 

8n*ann/350 . o 

ifCian eq.an) goto C70 

goto £80 

270 wr i t e ( 2 , 1 1 > nn 

IT format ( 5x nn= ' 15 ,/ ? 

ak* .56’5 

wr It e ( 2 , 22 J ( t inan>»'na’»'ai/aK-2T3 ) 

£2 format ( 5x /' lower convective zone temperature =‘fG.c,/3 

£8 0 gradi J * ( T xnd ( 2 )“T xnd( 1 ) )/del xnd 

9 rad 2 2= ( t xnd( i )“t xnd( i-1 ) ) /del xnd 
grad 2 3= C t gnd ( 2 ) -t gnd ( 1 ) ) /del xnd 
I 1 1 cont inue 

close( unit=l 3 

St op 

end 
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C INrUT DATA TILE C 

c roE c 

C ETFECT OF VARI0U5 fARAflETERE ON THE rEErOrrnANCE C 
C OF SALINITY 9EADIENT SOLAR FOND C 
C C 

c c 
c AFP 2 MDIX D TOSESH K. curTA c 
c ,- 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

t .5 0 .E 1 . AL> AL1 , ALE 
35.6AI .569 .519 HC,AK,AKS 

I.A61e-0T 1.035T®-06 EES ALPHA, ALPHAG,HA 
0.95 AE51, 

7.0 EC . 5 PS , .-^Hr 
C I 
7 0 J 

6T60 0.1 3600 365 NA , OELX , DELT , I CP 

366 366 NWE,NHA 

AVC SOLAR RAD. IN JAN , FEB , HAR , APR , MAY , JUN , JUL , AU6 , SEP , OCT , NOV , DEC 


0 . 

4 . 96 E71 . 433 . 556 678 . 

615. 55S. 430. 268. 99 

. 5. 0. 

0 . 

EO . 763. 364. 

540 . 

670. 7E7. 

7E8. 659. 

534. 358. 

765. E7 . 0 . 

1 . 

66. S60. 475. 

655 . 

787 . 84 

5 

844. 769, 

635. 458. 

E49 . 64 . 7 . 

7 0 

. 730. 339. 554 

. 7S7 

548 . 

97 0 

. 903. 833 

. 699. 5E3 

378 ICE 70 

30 

, 777. 3C3. 5C7 

. 746 

865. 

9E5 

. 9E0. 644 

. 779 547 

35E 7 58 E7 

36 

, 764, 339. 510 

657 

. 758. 

809 

. 873. 750 

. 639 483 

378 759 36 

S6 , 

, 730. E75 . 4S5 

, 549 

. 640 . 

67E 

. 677. 606 

. 578 47E 

E73 747 3£ 

7 E . 

707. £43 , 39E 

57 E 

. 607 . 

640 

. 659 . 595 

. 510. 39E 

E55 707 73 

3 . 

79 . E4C . 431 . 

587 . 

698, 75 

7 . 

747 . 689 . 

577 . 477 . 

E44 74 3. 

0 . 

36. EOE 406. 

55 7 . 

701 . 764 . 

758. 69E. 

569 . 391 . 

791. 33 0 


0 . 9 . 

1 28 

. 31 9 

495. 

620 

. 685 

680 . 

613 

. 490 

. 315 

1 25 

10 0 . 


0 . 3 . 

85 . 

259 . 

4E5 

543 . 

605 . 

605. 

538 . 

41 8 . 

255 . 

86 . 

3 . 0 . 


AVE APIS T 

EHP , 

IN JAN 

: , FEB 

; MAR , 

APR, MAT, JUN, JULj 

, AUG, SEP, OCT, NOV, DEC 

10.9 

10.4 

1 0 0 

9 6 9 

. 2 8 

.9 8. 

6 8 7 

11.0 

13.5 

16.3 

17.0 

17.6 

19.6 

7 9.7 

7 9.5 

1 9 1 

1 T . 2 

15.5 

14.1 

13.2 

12.6 

11.9 

11.3 




14.6 

14.1 

1 3 5 

13,0 

12 . T 

12.3 

12.0 

12.4 

15.2 

IT . 7 

17.2 

21.7 

23 . 0 

E3 . 6 

23 . 9 

23 8 

23 . 4 

EE . 0 

20 . 0 

18.6 

17.5 

16.8 

16.0 

15.4 



12 5 

18.9 

18 4 

1 T 7 

13 3 

16.8 

16.6 

18.0 

£1 . 1 

23 . 6 

25 . 7 

27 . 1 

28 . 0 

E8 . 6 

28 9 

28 . T 

2 8 , 4 

27 . 2 

25 . 3 

£3 . 7 

22 . 6 

21 . 5 

21 . 0 

17.3 



E5 . 6 

24 . 8 

E4 2 

2 3.5 

22 . 9 

22 . 6 

22 . 6 

24 . 9 

28 . 1 

30 . 5 

32 . 2 

33 . 5 

34 . 3 

34 . 8 

35 . 1 

34 . 9 

34 . 7 

33 T 

31.9 

30 . 2 

29 . 0 

28 . 0 

ET . 4 

26 . 5 



E9 . 1 

28 . 4 

27 . 8 

2T . 2 

26 . 7 

26 . 3 

2T . 1 

27 , 3 

31.8 

34 . 1 

35 . T 

36 . 9 

37 . T 

38 . 3 

38 . 5 

38 7 

38 8 

37 . 5 

35 . 9 

34 . 1 

32 . 6 

31 . 5 

30 . T 

27 . 0 



31 . 7 

31.2 

30 7 

30 4 

30 . 0 

29 . 6 

30 . 0 

31 .4 

33 . 6 

34 . 7 

35 . 9 

36 . 8 

37 . T 

35 . 3 

38 . 5 

38 . 4 

38 . 2 

37 . 7 

36 7 

35 . 4 

34 . 3 

33 . 6 

32. 9 

32 . 3 



E9 . 7 

29 . 4 

29 , 1 

28 . 8 

28 . 6 

28 . 4 

28 . 6 

29 . 3 

30 . 5 

31 . 2 

32 . 0 

32 . 8 

33 . 3 

33 . 6 

33 . 7 

33 7 

33 . 4 

33 . 0 

3£ . 4 

31 .5 

31 . 0 

30 . 5 

30. 1 

29 . 9 



27 . 9 

27 . 7 

27 . 5 

27 . 3 

27 . 2 

27 . 1 

27.2 

28 . 0 

26 . 5 

29 . 5 

30 . 3 

30 . 9 

31 . 3 

31 . 5 

31.6 

31 5 

31.4 

37 . 0 

30 .3 

£9 . 6 

£9 . 1 

28 . 5 

£8.4 

£8 . 1 



£6.7 

26.4 

26 . 1 

£5 . 9 

£5.7 

25.5 

£5 . 6 

26 . 5 

28 . 1 

29 .E 

30.3 

31 . 1 

31 . 6 

31 . 9 

32 . 0 

31 . 8 

31.5 

30 . 7 

£9.5 

£8 . 7 

28.2 

£7.7 

£7.3 

£6 . 9 



EE. 6 

EE . 7 

218 

£1.5 

£1 . 1 

£0.8 

EO . 8 

22.3 

25.2 

27.0 

£8.5 

£9.8 

30 . 4 

31 .£ 

31 . 3 

311 

30 . 6 

£9 . 0 

27. E 

E5.9 

25. 0 

24 . 4 

23.5 

23. 0 



16.5 

16.0 

15.5 

15.1 

14.8 

14.5 

14.3 

15.3 

18.9 

21 . 6 

23.7 

£5.£ 

25.9 

E6,£ 

£6.4 

26 . 0 

E5 .S 

EE. 9 

E0.9 

19.6 

16.7 

18.0 

17.3 

16.8 



7 1.5 

110 

10 6 

10.3 

10.1 

9.7' 

9.5 9 . 

8 12 

.3 14. 

9 17. 

3 19 

. 0 20 . 

I 

EO . 7 

EO . 9 

20 7 

19.7 

1 7 . E 

15.9 

14.7 

13.9 

13 . 2 

12.5 

11.9 



eO 81 

83 

85 83 

86 87 

87 ■ 

79 67 

57 49 

44 41 38 

38 40 

47 ! 

57 65 

69 72 

75 77 

AVE 

. RH 

IN JAN 










65 68 

70 ' 

72 73 

75 76 

73 < 

63 53 

44 38 

34 31 29 

£8 £9 

33 ■ 

41 48 

53 57 

54 6E 

AVE 

. RH 

IN FEB 










56 59 

61 

64 66 

66 59 

66 ( 

64 46 

38 34 

30 £7 £6 

£5 £5 

£5 : 

33 39 

44 47 

50 53 

AVE 

. RH 

IN MARCH 









40 4E 

45 ■ 

47 48 

5 1 52 

47 39 34 

30 26 

23 2 1 19 

19 19 

£0 i 

23 26 

30 3£ 

35 37 

AVE 

. RH 

IN APRIL 










37 

39 

41 4£ 44 46 

4 6 

43 

37 

34 

30 

26 

24 

22 

£0 

20 

1 9 

19 

22 

25 

29 

32 

33 

36 

AVE. RH I'N 

MAY 
















47 

49 

51 53 54 56 

56 

54 

48 

45 

41 

37 

34 

3E 

31 

30 

30 

30 

32 

33 

39 

41 

44 

46 

AVE RH IN 

JUNE 
















7 5 

78 

79 80 81 8£ 

61 

T9 

73 

68 

64 

61 

59 

56 

56 

57 

57 

58 

61 

66 

69 

72 

74 

75 

AVE. RH IN 

JULY 















87 

88 

89 89 90 90 

9 0 

88 

80 

77 

73 

7 0 

66 

67 

66 

67 

67 

69 

79 

78 

61 

83 

85 

86 

AVE RH IN 

AUG 
















ft3 

44 

85 85 86 86 

87 

64 

74 

66 

63 

58 

56 

54 

53 

54 

55 

59 

66 

7! 

75 

78 

80 

81 

AVE, RH IN 

5EP1 
















73 

75 

77 78 79 80 

80 

76 

63 

54 

48 

42 

36 

36 

34 

34 

37 

43 

51 

56 

63 

66 

69 

70 

AVE. RH IN 

OCT . 
















70 

73 

75 76 78 79 

60 

78 

64 

53 

43 

39 

35 

32 

31 

31 

34 

42 

SO 

56 

60 

63 

66 

67 

AVE RH IN 

NOV. 
















"79 

80 

81 63 64 6'5 

66 

66 

76 

65 

56 

49 

43 

40 

36 

38 

41 

50 

60 

66 

70 

73 


'75 77 AVE . 'RH IN, DEC 
0. LOAD 



